UMTG - 2 



Supersymmetric Biorthogonal Quantum Systems 

Thomas Curtright^'"'"'*, Luca Mezincescu^'*, and David Schuster^'* 

Department of Physics, University of Miami, Coral Gables, Florida 33124 
' ' School of Natural Sciences, Institute for Advanced Study, Princeton, New Jersey 08540 



February 1, 2008 



(N 

w): 

<. 

Abstract 

^ ; 

^ We discuss supersymmetric biorthogonal systems, with emphasis given to the periodic solutions that 

occur at spectral singularities of PT symmetric models. For these periodic solutions, the dual functions 
fSJ ' are associated polynomials that obey inhomogeneous equations. We construct in detail some explicit ex- 

^ , amples for the supersymmetric pairs of potentials V± (z) = —17 (z)'^ iLz-j^U (2) where U (z) = X]fc>o ''^kz^ . 

■ In particular, we consider the cases generated by U {z) = z and z/ (1 — z). We also briefly consider the 

effects of magnetic vector potentials on the partition functions of these systems. 
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1 Introduction 



' There has been some recent theoretical interest in non-hermitian Schrodinger equations, in the guise of "PT 



symmetric theories" [T]. Several previous authors have considered supersymmetry in this context 2J. Here, 
we consider a few elementary soluble examples and explore them in some detail. We believe exact solvability 
permits the underlying structure to be appreciated more completely. Some recent papers [51 [7] have also 
touched on the relevance of biorthogonal systems 4^ for PT symmetric models. As in [7], we wish again 
to stress the importance of such systems, and their generality, only here in the context of supersymmetric 
examples. We are not aware of any previous systematic discussion of supersymmetric biorthogonal systems 
along the lines of that given here. 

We consider models whose Hamiltonians are of the form 

C/(z) = ^t;fez'= (2) 

fc>0 

where the exponents k in the "superpotential" U are of the same sign, and where v and Vk have arbitrary 
values. We look for eigenfunctions and associated functions that have particular analytic behavior near 
z = 0. The connection to PT symmetric theories is achieved by writing z = me*^, as explained in jJJ- 
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2 General Theory for u = 

First we consider — 0. For bounded functions of a: G (—00, +00), with z = me", the spectrum of any 
Hamihonian of the form O [7] 

^ ^ fe>0 

is known to be the real, positive half-hne, for any choice of fis such that X]fe>o l/^'^l ^ ^-"^^ '^'-'^ J^^^ ^'-"^ 
those ^s which are real or x-translationally equivalent to real values. So PT symmetry is not required for 
real energy eigenvalues in these models [S]. Here, we will restrict our attention to 27r-periodic functions 
of X, with z = me", and their duals, to obtain a discrete subset of real energy eigenvalues, namely just 
{En^n^ I n = 0,l,---}. 

In this situation the general theory for supersymmetric pairs of non-hermitian Hamiltonians goes as 
follows. For a given U (z) we may construct pairs of finite polynomials in z~^, 

n 

x^W--;rE<.-^' (4) 



that satisfy the inhomogeneous equations. 



z^^±U [z) ) {z) + nxt (z) = A± {z) ^ J] A± (k) 7? (5) 

fe>0 



where U {z) and (z) are analytic about the origin, at which point they all vanish. A priori the A^ (z) 
need not be given, for a given U (z) they may be determined along with xt, (-^)- Clearly, xt (^) c^,o/-^" 
and, by convention, we normalize so that the coefficient of the most negative power of z is just c^g = 1. 
This choice also fixes the normalization of A^. For a given U [z) the four functions Xn {z) and A^ (z) are 
now completely specified. 

Alternatively, we construct pairs of functions i/)^ (z) analytic about the origin, hence given by series in 
non- negative integer powers of z, that satisfy the equations 

z^±C/(z))^±(z) = n^J (z) (6) 

with 7/;^ (z) z". Usually these are infinite series, but again, for a given U (z), both fimctions V'n {z) 
are now completely specified. 

It follows that the (z) are eigenf unctions of H± with eigenvalues n^, 

H±i;t (z) = (z) (7) 

d , A / d , , A / d \^ , 2 . d 



H,^^z-TUiz)j^z-±Uiz)j^^z-j -Uizr±z-Uiz) (8) 

while the associated polynomials obey the inhomogeneous equations 

{H± k') (z) - (^^ T C/ (z)) A± (z) - kA^ (z) (9) 

Note the RHS of this last equation involves only positive powers, z" for n > 1. Moreover, it also follows 
that these functions are biorthonormal systems such that 

±-J^xUz)i^tiz)^S,„, f;;^±HV,±(z)=^ (10) 

n=0 

These relations encode orthonormality and completeness on analytic functions. For more or less obvious 
reasons, we will call the full set of fimctions {xj' j'4'k:Xr ^'^n} analytic, supersymmetric biorthogonal 
system corresponding to the given function U (z). 
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The dual polynomials and their inhomogeneities. We start the construction with the dual polyno- 
mial^ and explicitly compute the first few. This serves to illustrate how the various functions are uniquely 
determined, given that xt i^) We choose the coefficients in xt to eliminate all z^*^ terms on 

the LHS of ([5]), for fc = 0, 1, • • • , n, and then we sweep all the remaining positive powers of z into the (z). 
We find 

X and A Table 

Xo^W = l Ao±(z)=±C/(z) = ±Efc>o^fcz' 



xt (^) 



V 



z z I A fz^ 'V^ I 



etc. When the sum X]fe>o ^^^^ defines C/ is finite, the process is clearly finite mathematics all the way, for 
any n. Note that Xn (-^) ^'^d A~ (z) are obtained from x+ (z) and A+ (z), or vice versa, just by flipping 
the signs of all the vs. Also note that all negative powers of z can be expressed as finite sums of either 
{X+ (^) I ^ > 0} or {x- (z) I n > 0}. 

More systematically, we solve ([5]) as follows. We impose the condition that the RHS involve only 
positive powers of z so as to obtain an inhomogeneity that will be orthogonal to the span of {z" | n > 0} 
under contour integration §^ ■ ■ ■ . This leads to n equations that fix the coefficients c„^fc for /c = 1, • • • , n 
in terms of Cnfi, the latter being an overall choice of normalization. Thus 

n n n 

^ c±^. (j - n) z^ ± ^ Vkz'^ Yl + " E = E ^" (^) (11) 

j=0 k>0 j=0 j=0 k>0 

j 

SO (i -n)±Y^ VkC^^j-k + =0 for j = 0, • • • , n (12) 

k=l 

n 

and Xt (k) = ±J2^k+n-jci^, (13) 
We determine all the coefficients in the xt i^) from the pair of equations 

where c^n = 1- Now det ( . " 1=7 (2n - ?) > for < 7 < n, and I ^ ^ . 

"''^ \ n j-n J ^ ■" j-j \ 



J ~ n n 



j{2n-j) \ _^ n - j 



So we have the recursion relations 



±1 x-^ 



Each ct j depends on only the first j coefficients in the expansion for U (z), i.e. just on Vk<j. From the j 
we then determine the (k) using (fT3|l . Note that A± (fc) 7^ for /c > n is possible here, depending on the 



^So these systems are exceptions to the statement: "A direct and frontal attack on the problem of determining [the dual 
polynomials] is not usually fruitful." - Morse and Feshbach [4] p 931. 
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values of the vs. The (fc) will depend on all the vs, in general, with k taking on all values up to and 
including the highest power of z appearing in U. 
For example, for n — 2: 



^2,0 



-2,1 



-2,2 



(k) = Vk+2 + Vk+lC^j + VkC2 2 



Vk+2 ± -^ViVk+l + ( ±2^^ 



±2-2 - -^v. 



1 



6 



(15) 
(16) 



And for n 



-3,0 



± , 1 1 2 

C3,2 = ±41^2 - Yo«l 



ibAg (/c) = + Ufe+2C3_i + Wfc+lC3^2 + -feC3_3 



± 3f lW/c+2 

5 



- ) vu+i 



3,3 33 20 ^ 30 ^ 



±\v, - ^^v,v, T ^vf ) Vk 



(17) 
(18) 



Hence the table given earlier. 



Energy eigenfunctions The energy eigenfunctions, when written as series. 



(19) 



can be determined just as the dual polynomials were by direct solution of ([6]), or else the eigenfunctions 
can be determined by imposing the bi-orthonormalizations in (|10p . As discussed in [7], these orthogonality 
conditions amount to a set of triangular equations which can always be solved, sequentially, for the j in 
terms of the c^s. Namely 



1 



k — n 

5]4fc-„-j««,j =0 for k>n 



(20) 



The series for is a development in the minors that invert these triangular equations. By considering all 
k > n in succession, we obtain all ^ in terms of c^;, or vice versa. 

For convenience, we again choose the normalizations c,^q — 1. The results of the recursion relations ([20|) 
are the pair of correlated series 



1 



1 

'^n+2,2 
'^n+3,3 



Z 
1 

''n+2,1 
'-n+3,2 




1 

^n+3,1 



1 "^n+l,!^ + 



c: 



z 


1 



-n+2,2 '^n+2,1 



C± C± C± 

^n+3,3 ^-71+3, 2 '-n+3,1 



z^ + - 



(21) 



(22) 
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Alternatively the c s may be expressed in terms of the a s. 



(^) 



Xn [z) 



-*n-l,l 



^n,2 
Z 
1 



a. 



± 



n-2,2 



a. 



± 

n — 3,3 



a 



''l,Tl-l 



— 1 - a^ 



n-3,2 
1^ + 



a 



a 



+ -••• + (-1)" 



^n-2,2 



''n-3,3 



n,3-= 



1 

^n-3,1 



l^n-3 
± 

n-1.1 
± ■ 
n-2,2 

1 



Z 




1 



'^0,n-3 
1 

^n-2,1 



a. 



n-2,1 



z'- 



a. 



n-3,2 





1 

*r!-3,l 



-'n-2,2 
^ri-3,3 







1 



n-2,1 
^n-3,2 




1 

^n-3,1 



(23) 



(24) 



Orthogonality and Completeness With either pair of these correlated series, the orthogonality relations 
in pI7|) are easily checked. They amount to the obvious statements that 



1 



dz 



Xn (z) iz) = l , 



1 



dz 



X„ (^)C+fe(^) = for fc>0 



(25) 



as well as more involved cancellations to show 5^/^ Xn+fc (-^) C-^) = for A; > 0. In particular, 
there is a complete cancellation of all the contributions to this latter contour integral upon expansion of the 
{k + 1) X {k + 1) determinant that appears in the first line of p4|) . Exploiting the multi-linearity of the 
determinant, and performing the integrations entry by entry. 



1 
27ri 



dz 



xUk WV-^ {z) 



a 



1 

z 
1 



dz 



a 



n+fc-2,2 '^n+fe-2,1 
^n+fc-3,3 '^n+fc-3,2 



'^0,ri+A; '^0,n+fc-l 



a 



z 

1 

n+fc-3,1 



^l^n+fc-3 
\,n+k-2 






1 



-'0,n+fc-3 



^ + an,2^ +a„,3^ 



^n+/c— 1 







-^0,1 









^n+fe-2,2 



1 

'^n+fc-2,1 



^n-\,k+\ ^n-l,k 



a 



i,k-2 



1 

n+fe-3,1 



"'n,fc-2 

a'n-i,k-i a 



3 



1 



"'n,fc-3 
± 

n-l,fc-; 



a«,o = 1 






(26) 
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which vanishes since the 1st and the (fc + l)th rows are identical. Similarly 



1 



dz 
z 



^n+fc,fc-2 



= 



(27) 



for fc > 0. 

The correlations between coefficients in (pij) and (1^^ . or in ([^51) and ([M]) . also imply the complete- 
ness relation in (jlOp by guaranteeing that all terms of the form z^w^^ for fc 7^ / cancel out in the sum 
SJ^o (^) {^)- This cancellation is encoded in the identities 



+ E E + E E E <.<.-"^'-" 



(28) 



n=0 i=l 



n=l k=l 



n=l 1=1 k=l 



which follow from (I20p . Terms of the form {z/w)^ as required to give jirx ^'I'e provided just by the leading 
terms in ^ and ([22l) . 

in ^/;^ (z) are again finite polynomials in the vs. While convergence of this 



The coefficients of 



yU+k 



series, as written, is certainly not obvious for arbitrary ws, it is clear that convergence can be determined on 
a case-by-case basis from the explicit form of the coefficients. Moreover, when the number of vs is finite, 
no matter what their values are, it is not too difficult to show the tp'^ (z) are entire functions of z. Thus, 
if we assume the requisite convergence of the series, all the {-0^ (z) \ n> 0} are determined, and either 
(z) I n > 0} or (z) | n > 0} are complete on the span of {z" | n > 0}. All positive powers of z 
can be expressed as series of {"0^ (-2^)}? just as all negative powers of z can be expressed as finite sums of 
{Xt {z)}- 

Remarkably, the non-degenerate energy eigenfunctions {ip^ (z)} just obtained turn out to be all of the 
eigenf unctions of which are 27r-periodic in x, where z = me*^. Moreover, the fact that (z) are 
indeed eigenfunctions, as given in ([7]), can be deduced in a novel way from the biorthogonality of 
{xf {z) li'k (-^)}' ^'^^ completeness of (z)} for analytic functions about the origin, as described in 
[7]. In fact, the argument given in [7] can be adapted to the ffi'st-order equations. Completeness on analytic 
functions about z = allows us to write 



z-±U{z) 
dz 



00 

E 

k—n 



(z). 



(29) 



Note the chosen interchange V'^ ^ 4'^ upon LHS^RHS. From this expansion and biorthonormality, we have 
= ^f^Xk (•^) {zj^ ^ ^ i^)) V'n (z)- But then, upon integrating by parts and using ([5|) as well as the 
orthonormahty relations, we also have -^f^xt i^) ^ ^ (^)) V'n (z) — n6k,n- So hn,k = nSk,n, and ([6]) 



is obtained!! Conversely, given 
for fc^ ^ ^ just by inserting 



and (O, we may prove the orthogonality relations §^Xk i^) ipn (z) — 
± [/ (z)) and integrating by parts. That is to say 



n^^xi (z) i>t iz) = §^xi {z) [z^^ T U (z)) (z) 



z 



-^^tc/(^))xJ {z) 



(30) 



Thus {n^-^e)§^xi {z)^tAz) = Q- 

^Note that 2 = is a regular singular point of Jg} and/or J?}, for any number and any values of the vs such that J^j-^q ^kZ^ 
converges near the origin. In fact, II22I I is exactly the conventional series obtained by expanding about the regular singular 
point at 2 = 0, albeit the series was obtained here in an unusual way from the properties of the dual space polynomials. 
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3 Examples 

As an explicit example, to parallel the discussion in [7], we note that the superpotentiaH 



gives a simple quadratic potential for 



V+ = + z-^U = (32) 
az 



but a much more complicated partner potential for iJ_ 



The complexity of iJ_ suggests that we seek a simpler J7 to fully illustrate the general theory. 
Complex Morse potentials Again referring to [7], we consider U (z) = /iz, hence 

H± = (^z-^ + ± fiz - ^^z^ (35) 

Note that for this simple example the solution of one Hamiltonian, say immediately gives the solution 
for the other, through the relations [z] — [—z], ip~ (z) oc "ip^ {—z). But this is not necessarily the 
most transparent way to write the solutions for H_ . 

When the vector potential is not present, v — it may be best to simply note 

^ ^) (^"-1/2 (^) =F In+l/2 (z))^ = n (/„-l/2 {z) ± In+l/2 {z))^ (36) 

So then it is obvious that 



are eigenf unctions of 



iz) = Zt\j- {ln-l/2 (Z) T 4+1/2 iz)) (37) 



d \ ^ o f d \ f d 



^-n^^j -Vd-z^')['d-z^v ^''^ 

with eigenvalues rt^ as given by 

H±^^ (z) = nV^.t (z) (39) 
and with normalization constants Z^. Other ways to write the eigenfunctions for this example are: 

^ii^)-^J4Tz + l + n)'-I^±^ (40) 



V2 \ dz 

*i (=) = (1)" it WTtTTTn^ (I)" ^£ SFTTTTT^ (1)"") 



If we choose Z+ — then the relations between the two sets of eigenfunctions are just 

dz 



z-^ ± z ) V± (z) = n^J (^) (42) 



''Actually, for any constant c the supcrpotential U {z) = z y. j^^^j^^y^|^j gives a simple quadratic potential for Hj^ 
= —U^ + z-T^U. But if c ^ 0, (7 involves a logarithm. In 2, and hence is not a periodic function of x for z = me'^. 
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Now, what about the dual polynomials {xt \ n> 0} which are the biorthonormalized duals for \ n > 0}? 
These are given by 

^ ^ (!)■■ ( E ^^^^^^^ . ^^^n^ <«) 

As illustration of the general theory, for all fc, n > we have the orthonormality and the completeness 
relation as given in (fTO|) . It is straightforward to check these relations by using the explicit series forms for 
xt and ip^. 

Again choosing both normalization factors to be the same, = Z~ = Z„, the dual polynomials for the 
complex Morse potential are solutions to the exceptionally simple inhomogeneous pair of equations 

± ^ j xt (z) + nxt {z) = zX^ (44) 



where by direct calculation we find 

~ z„ r(K2j + i) 



.J.(4!^!!LM|±il, .J (45) 



That is to say 

= ±<„ (46) 

Moreover 

{H± - n") xt, (z) = (^z-^ T z^ zA± - znX^ = z (A± - nAj) T ^'A± (47) 
or with the explicit coefficients 



f_n"+^+L"/2J r (In/21 + i) 
(H- - n") ..T ^ '"^^ (d + (-!)"») ^ + (49) 

The coefficients are a bit awkward, particularly the phases, but are dictated by xt (z) ~o 

Singular potentials Now we go on to discuss models with several /iS. In particular, if the sums over k 
are infinite, the potentials can have fixed singularities for finite values of z. We explore the situation for a 
particular supersymmetric pair of such singular potentials. Namely those generated by the superpotential 

[/(z)==^ (50) 
Up to the scale of z, this is the unique superpotential that reproduces itself to obtain V+ (z) = U (z). 

What are the exact energy eigenfunctions for these potentials, analytic about z = 0? First consider the 
Hamiltonian with potential V+. 

H 2 '^^ ^ z ^ \ ( + ^ \ (52) 

dz"^ dz I — z \ dz 1 — z J \ dz 1 — z/ 

(z) = 2;" (1 - z) hypergeom ^ 1 + n + \/ \ + -n? , l + n- \/ \ + -n? , [1 + 2n] , z^ (54) 

{z) = l-z 
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We note the especially simple form for the ground state. Excited states are not such elementary functions. 
The dual polynomials in this case are 



\ fe=i 



(2n-A:-l)! r{k-n + V^?TT) T {n + Vn^ + l) 



fc! (2n-l)! r (1 - n + v^^^TT) r (-fc + 1 + n + Vn^ + l) 



xtii) 



r (1 + 2n) r (i - n + VTT^) 

These are solutions of the inhomogeneous equations 

{H - n^) xt {z) = V+ (z) x+ (1) = X+ (1) 



(55) 
(56) 

(57) 



where the coefficient of the singular inhomogeneity is just xt, evaluated at the singularity. The orthonor- 
mality relations between the eigenfunctions and the dual polynomials are again the expected ones, (|10p . In 
this case, the contour encloses the origin once in the positive counterclockwise sense, but lies within the 
unit-radius circle of convergence of the series for -0+ [z). (Or at least, the contour swerves "to the left" to 
avoid the singularity at z = 1.) From the explicit series there also follows the expression for the Cauchy 
kernel as given in PH)) . 

Now consider the Hamiltonian with the superpartner potential, . 



dz'^ ~'~ dz 



z + l 
{l-zf 



d z 
z-r + 



2 



dz 1 — z 
z + l 



dz 1 — z 



dz 



dz 



(l-z) 



2 V-n =n 



ifj^ (z) = z" (1 - zf hypergeom 2 + n + \/l + rfi, 2 + n- \/\ + "n? , [1 + 2n] , z^ 



(58) 

(59) 
(60) 



1 - z 

Again, we note the especially simple form for the ground state. And again, excited states are not such 
elementary functions. The dual polynomials are now given by 



Xq (z) = 1 and for n > 



(61) 



After simplification of this last expression, we obtain 



1 - z 



(x+W-x+(i)) 



1 



fe=i 



z 
"fc! 



n + 2nfc — fc 



(2n-fc-l)! r (fc - n + Vn^ + 1) r (n + + l) 



(1 - z) , 2 in - j) r (2n - j) T (l + j ~ n + V^^^^Tl) T (n + \A?TT) 

— y z-' - 



(2n-l)! T{l-n + Vn'^ + l) T (-fc + 1 + n + Vn^ + l) 

(62) 



j!(2n)! r (1 - n + %A?TT) r (-j + n + Vn^TT) 



Note that Xn (1) — 0- J^^^ as ([50]) has an additional factor of (1 — z) compared to ([5l)l . so too does (IS^ 
compared to ([55)) . 

These results also exhibit biorthonormality and lead to yet another expression for the Cauchy kernel, 
as in (jlOp . as may be established from the explicit series. The inhomogeneous equation obeyed by Xq is 
obviously just 

H-Xo - V^Xo - (63) 
(1 - z) 

On the other hand, the inhomogeneous equations obeyed by Xn>o ^ more interesting. 

z 



{H-n^)x- {z) = V+{z)2x-'{l)^- 



d 



(64) 
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This is the same singular inhomogeneity as appears in the superpartner dual polynomial equation, ([57]) . only 
with a different coefficient. In fact, 

2^X^(1) = -nX^(l) (65) 

All this is more transparent upon exploiting the Darboux factorizations of the Hamiltonians, (|52p and 
As expected from the general theory, we have 



(i/± - n^) V± (z) = (66) 

- n^) xl (z) = Vo [z) x+ (1) = (1) (67) 

- n^) x~ (z) = -l/o (z) nx+ (1) = Y^^'^^n (1) (68) 
bearing in mind (j56p and (|65p . as well as 

zj-^±^^^^{z) = n^^{z) (69) 

+ ) (^) = (^) + 7^*^"° 

The last of these is strikingly simpler than expected, exhibiting an inhomogeneity only for the case n = 0. 
Then there are the complementary, but subsidiary, relations: 

which allow the functions rj^ (z) to be constructed from ([5^ and ([55]) . Differentiating again gives 

(i?^ - n^) {z) = -nv^ (z) - ^~zj- T ^) 77± (z) (73) 

4 Magnetic field effects 

Now consider v ^Q. For a given [/ (z) = X]*:>o ''^fe-^'' '^^ define 

//, = (4 + .) ' - 1, ± = (4 + ^ T t/) (4 + . ± c/) (74) 

Obviously, H± [v] = H± [— i^]. We seek eigenfunctions and associated functions that have particular analytic 
behavior near z = 0. 

Energy eigenfunctions We look for pairs of functions (z) that satisfy the equations 

z^ + i^±f/(z)^V'^(z) = (n + i.)7Aj(z) (76) 

with (z) z". Usually these are infinite series, but again, for a given U (z), both functions (z) 
are now completely specified. We note that negative integer n are now admissible, and independent of 
the corresponding for generic i>. It follows that the (z) are eigenfunctions of H± with eigenvalues 

iJ±^±(z) = (n + i.)'^^ (z) (77) 
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It also follows for generic v that a suitable dual to ■0^ (z) is provided by irl^ (z), where H±iIj^j^ = 
{—n — 'i/'^n so that ip^j^ has the same energy as ip^. Nevertheless, as we discuss shortly (cf. remarks 
about the "right-sector" given below) , if we consider subsectors of the spectrum, we may use polynomials as 
alternative dual functions, just as in the periodic situations discussed above for i> ^ 0. 
For example, reconsider the singular potential U — For generic v the solutions of 

which are single-valued about z — are 



(z) = z" (1 — z) hypergeom 



[l + 2n + 2v] , z 



(78) 

where the indices on the hypergeometric function involve p, the roots of = —1 — — 2nv — 2vp + p'^. That 
is to say p ~ v ± \J 1 + {n + v)^. The other solutions of the equation are 

^-n-2i/ ^-j^ _ hypergeom {[-p + 1 - + p ~ n - 2v\ ,[\ - 2n - 2v\ , z) 

but these have a branch point at z = and are not single-valued unless 2v e Z. 
For the partner potential, the single-valued solutions of 

are 

-0- (z) = z" (1 - zf hypergeom ^ 2 + n + z/ + ^ I + {n + vf, 2 + n + v- \J I [n ^ vf ,[l + 2n + 2v] , z^ 

(79) 

while the solutions with a cut, for generic i/, are 



Z 



(1 — zf' hypergeom {[2 + p — n — 2v, —p + 2 — n] ,[1 — 2n — 2v] , z) 



where p = v±^\-\-[n + vf as before. 

The dual polynomials and their inhomogeneities. On the "right sector" {i/i^ (z) \ n> 0} (this ter- 
minology is explained in [7]) we also look for pairs of finite dual polynomials in z^^, of the form given in ([4]). 
These satisfy the inhomogeneous equations 



z^-v±U{z)] xt {z) + [n + I.) x,T (z) = A± (z) ^ ^ A± (k) z^ (80) 

fe>0 



dz 



where U (z) and (z) are all analytic about the origin, where they vanish. A priori the A^ (z) are not 
given, but are to be determined along with xt i^)- Once more we normalize these polynomials so that 
Xn {z) jr;Q 1/z". This choice also fixes the normalization of A^. For a given U [z) the four functions 
Xn {z) and A^ (z) are now completely specified. 

The associated polynomials also obey the inhomogeneous Hamiltonian equations 

H±xi {z) = (^^ - ^ T f/) - ± C/) xt (^) = (^^ - T C/) {^■t (2) + xt (^)) (81) 

= (^^ - ^ T C/) A± (z) -{n + v) {Kt (z) - (n + ^) xt (^)) (82) 



That is to say 



(h± -{n + i^f) xt (z) = (^^ - ^ T C/ (z)) At (z) -(n + iy) (z) (83) 
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Note the RHS of this last equation involves only positive powers, for n > 1. 

Moreover, it also follows that for k,n> and —l/2<i/<l/2 these functions are biorthonormal systems 
as in PU]). These relations once again encode orthonormality and completeness on functions analytic around 
z = 0. For example, it is straightforward to show {n + v)^ § ^ Xk (-^) (-^) — {k -\- v)^ § ^ ■0^ (z) xt (z) 
by inserting the Hamiltonian (17411 and integrating by parts to convert to ([75]) . That is 



= 



((fc + -{n + lyf) j> ^ xt (z) V',f (z) = {k - n) {k + n + 2i.) I ^ x^ {z) i>t (z) (84) 



dz , , X , , ^ ^ f dz 

z 

For k j^n and both > 0, and ^ Z<o, we conclude § ^ xt (-^) V'n {z) = 0. 

Once again the various functions are uniquely determined, given that Xn {z) j^;^ 1/z". We choose 
the coefficients in Xn {z) to eliminate all z~^ terms in (|80p . for A; = 0, 1, • • • , n, and then we sweep all the 
remaining positive powers of z into the {z). 

z^-i^±U (z)^ xt (z) + {n + u) x^ {z) = At (z) (85) 

n 

U{z)^Y.^,z\ xtiz)^-Y.<,'' ^ A±(z)^^A±(fc)z'= (86) 

k>0 j^O k>0 

n n n 

{j-n- v) z^ ± ^ v^z^ + + E = E (^) (87) 

i=0 fc>0 3=0 j=0 k>0 



So then 



<,0'-"-^)±E^fc<^-fe + ('^ + ^X, =0 for j=0,---,n (88) 

k=l 

n 

and (fc) = ± ^ Vk+n-jctj (89) 
i=o 

Therefore the recursion goes like this. First we determine all the coefficients in the Xn (z) from the pair of 
equations 

wherec±o^l. Now det f " + " " ^' " M = j (2n + 2. - j) and f " + " " ^' " 



1 /J^^i^^ n + z^ 
—n — V 1^ + n — j 



, so we have the recursion relations 



J {2n + 2v- ]) ^ 



Each c„ J depends on only the first j coefficients in the expansion for U (z), i.e. just on Vk<j- Then from 
the j we determine the (fc) using ([89|) . Note that A^ (fc) 7^ for fc > n is possible here, depending on 
the values of the vs. For example, for n = 2: 

±A^ (fc) = Vk+2 + W/c+lC^i + UfccJ_2 

And for n = 3: 
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■^'^ 2i^ + 3 6 (j^ + 2) (2i/ + 5) 2 (2j/ + 3) (2i/ + 5) " 

±A^ (fc) = Ufe+s + Ufc+2C^i + ^^fc+lC^2 + "f^fec^s 

= "'^^^ ± 2^"^"^+^ + {2^^' ~ 2(27T5) ^0 

±1 (4z/ + 9) „, „, _ 1 „ 3 



2u + 3"^ - 6 (. + 2) (2. + 5) "^"^ ^ 2 (2. + 3) (2. + 5) "^ j "'^ ^^'^ 
Thus we have the lowest four dual polynomials and their inhomogcncities. 

(z) = 1 , Ao± (z) ^±U{z)^±J2 Vkz"" (96) 



fc>0 



{z) = -±^v,, Af (z) - 5] f ^^iz;,. ± w.+i) z'^ (97) 



k>0 

^ ,2 



^'^ = -^^ 2i7T3"^I ^ 2FTT)"^ W^f' ^^^^ 

(-) = E ( (^(^-^ T 2(2^7^3)^0 ^ ^ "^+0 ^''^ 

± / ^ 1 1 1/1 1 2^ 



z3 2iy + 5 z2 ^ V 2(!/ + 2) ^ 2(2i. + 5) 

^ 2^ri:3"^ ~ eFrir^TT^"^"^ ^ 2 (2. + 3) (2. + 5)"' 

(^2i:+3"3 T 6(.+2)(2.+5)"2Ul " 2(2.+3) (2.+5) ^1 J \ ^^^^^ 

+ (2(rT2)^2 T 2(2;^^?) + 2l^^^l"fc+2 ± J 

etc. Note that Xn i^) ^^'^ i^) obtained from Xn (^) ('^); or vice versa, just by flipping the 

signs of all the vs. These results reduce to those in the previous Table, upon setting — 0. 



The partition function is not analytic in ly For generic v 



oc 



= trace (e-^I^^l) ^ E 

n— — oo 

= e""^' (i9 M + [-Z/] - 1) (102) 

upon defining the theta functioijf] 

oo 

i?M = ^e-"'-2"^ (103) 

71=0 

whereas for v = 0, the partition function is 

oo 

Zo^^e-'"-^ =i)[0]^ 1.3863--- . (104) 

rs=0 

So then 

hm Z [v] = 2Zo - 1 7^ Zo (105) 

»0 

This would suggest a phase transition as an indicator for bulk systems governed by these dynamics. 

''Or in terms of Jacobi's functions, !?3 {z,q) = 5Z5^_cx) so i} [i^] + 'd [—i/] — 1 = (z = iu,q = 1/e). 
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5 Conclusions 



We discussed supersynimetric biorthogonal quantum systems along the lines of [7] , paying particular attention 
to the structure of non-hermitian systems with periodic solutions, for which cases the duals of the energy 
eigenfunctions are not simply related to the eigenfunctions by either complex conjugation or PT reflection. 
We worked out the general theory for single particle quantum systems, and we illustrated the general theory 
with several explicit exact examples. 

It remains to investigate many-body or field theoretic extensions of these supersymmetric systems, say 
by adapting the perturbative methods in ^8 on supersymmetric Liouville field theory, or by employing the 
powerful non-perturbative methods of conformal field theory [5]. This additional study is in progress, 
and represents one application of the formalism presented in this paper. There is a rich literature on 
Liouville and super-Liouville theory, models whose importance came to light in the work of Polyakov on 
string theory [10 , but which were subsequently developed much further in the context of conformal field 
theory and its applications to critical phenomena, as well as to subcritical string theory [11| . In particular, 
the super-Liouville correlation functions 12J have been shown to exhibit interesting analytic behavior in 
the exponential coupling constant, similar to the analytic structure of correlators for non-supersymmetric 
Liouville field theory [13j. The behavior of these correlators is related to properties of various WZNW 
models |14| . and there are particularly intriguing features that correspond to purely imaginary coupling 
constants - precisely the field theory extensions of the type of models discussed in this paper. 

As for other applications of supersymmetric biorthogonal quantum systems, say to non-relativistic situ- 
ations, an interesting possibility would be to consider driven/dissipative condensates as suggested in [7], but 
with additional fermions in the condensate |15j . This too is under study. 
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